clc
ime = '001_02_00_nhf_r1.wrl';

tic
A = ucitaj(ime);

N = 4000;
A1 = sazmi(A,N);

r = 20;
e = 0.15;

% racunamo point signature s reduciranim skupom tocaka
S1 = pntsg(0,10,-1790,A1,r,3*e);
% racunamo point signature s izvornim skupom tocaka, 
S = pntsg(0,10,-1790,A,r,e);

toc